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How unusual was autumn 2006 in Europe? 



G. J. van Oldenborgh 

KNMI, De Bilt, The Netherlands 



Abstract. The temperatures in large parts of Europe have 
been record high during the meteorological autumn of 2006. 
Compared to 1961-1990, the 2m temperature was more than 
three degrees Celsius above normal from the North side of 
the Alps to southern Norway. This made it by far the warmest 
autumn on record in the United Kingdom, Belgium, the 
Netherlands, Denmark, Germany and Switzerland, with the 
records in Central England going back to 1659, in the Nether- 
lands to 1706 and in Denmark to 1768. The deviations were 
so large that under the obviously false assumption that the 
climate does not change, the observed temperatures for 2006 
would occur with a probability of less than once every 10000 
years in a large part of Europe, given the distribution defined 
by the temperatures in the autumn 1901-2005. 

A better description of the temperature distribution is to 
assume that the mean changes proportional to the global 
mean temperature, but the shape of the distribution remains 
the same. This includes to first order the effects of global 
warming. Even under this assumption the autumn tempera- 
tures were very unusual, with estimates of the return time of 
200 to 2000 years in this region. The lower bound of the 95% 
confidence interval is more than 100 to 300 years. 

Apart from global warming, linear effects of a southerly 
circulation are found to give the largest contributions, ex- 
plaining about half of the anomalies. SST anomalies in the 
North Sea were also important along the coast. 

Climate models that simulate the current atmospheric cir- 
culation well underestimate the observed mean rise in au- 
tumn temperatures. They do not simulate a change in the 
shape of the distribution that would increase the probability 
of warm events under global warming. This implies that the 
warm autumn 2006 either was a very rare coincidence, or 
the local temperature rise is much stronger than modelled, 
or non-linear physics that is missing from these models in- 



creases the probability of warm extremes. 



1 Introduction 

Meteorologically, the autumn of 2006 was an extraordinary 
season in Europe. In the Netherlands, the 2-meter temper- 
ature at De Bilt averaged over September-November was 
1.6°C higher than the previous record since regular ob- 
servations began in the Netherlands in 1706, which was 
a tie between 1731 and 2005 (Fig. QJ). The excess is 
much larg er than the uncertainties in th e earlier part of 
the record: Ivan Engelen and Geurtsl (1 1985b estimate that the 
stan dard error of monthly tem peratures is 0.2-0.3°C (see 



also Ivan den Pool et all 1 19781) . The Central England Tern 



Correspondence to: 
knmi.nl) 



G. J. van Oldenborgh (oldenborgh® 



perature, 12.6°C, also was the highest since the beginning 
of the measurements in 1659, 0.8°C higher than the previ- 
ous record of 1730 and 1731. Pre-instrumental reconstruc- 
tions indicate that September-November 2006 very likely 
was t he warmest autumn sinc e 1500 in a large part of Eu- 
rope dLuterbacher et al.Ll2007l) . 

The impacts of the high temperatures on society and na- 
ture have not been very strong, as in autumn a higher tem- 
perature corresponds to a phase lag of the seasonal cycle. 
Flowering bulb farmers in the Netherlands were reported to 
have problems due to premature flowering. However, a simi- 
lar anomaly in summer would have given rise to a heat wave 
analogous to the summer of 2003, whic h caused severe prob- 
lems (e.g JSchar and Jendritzkyi 2004 ). 

In this article the heat anomaly of the autumn of 2006 in 
Europe is analysed. First the observations are shown and 
return times are computed under the obviously false assump- 
tion of a stationary climate. Next the first order effects of 
global warming are subtracted, and return times of the re- 
maining weather signal computed. The main weather factors 
are identified, and possible changes in their distribution are 
investigated using climate model simulations. 
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Fig. 1. Autumn 2-meter temperatures at De Bilt, the 
Net herlands with a 10-yr r u nning mean (green) 1706-2006 
(a, Ivan Engelen and Nellestiinl 1 19961) . 1901-2006 corrected for 
changes in th e thermometer hut, location and city effects (b, 
IVerbeedl2003h . 
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2 Observations 

In Fig.Q]two time series of autumn (September-November) 
averaged temperature in De Bilt, the Netherlands are shown. 
The first one covers 1706-2006 and is a combination of ob- 
servations at various loca tions in the Netherlands converte d 
to De Bilt temperatures (Ivan Engelen and N ellestiin. 
The second time series 1901-2006 has been corrected for the 
effects of changes in thermo meter hut, direct environment of 
the m easurement location (IBrandsma and van der Meulen , 
2007albh. and the effects of the growth of the cities in the 



area (Brandsm a et all 120031) . The value for 2006 is seen to 



be well outside the distribution defined by the other years in 
both series. 

Figure|2]shows two extrapolations of the cumulative distri- 
bution of the second, homogenised, time series 1901-2005, 
the first assuming a Gaussian (normal) distribution and the 
second using a peak-over-threshold method with a fit of the 
warmest 20 % to a Generalised Pareto Distribution (GPD; 
ColesL 120011) . These extrapolations are based on the obvi- 
ously false assumption that the autocorrelation is zero, i.e., 
that there are no long-term variations except those resulting 
from accumulations of interannual variability. From figure [2] 
it seems likely that the Gaussian distribution overestimates 
the probability of high temperatures. In spite of this a return 
time of 10000 year is obtained with this extrapolation. This 
shows that the assumption of no autocorrelation is very likely 
false: climate does change on longer time scales. Physically, 
global warming has made high temperature anomalies much 
more likely during recent years, and this and other long- 
term variations increase the probability of clustered high ex- 
tremes. 

The same analysis has been applied to all grid 
points of the 0.5° Global Historical Climate Network 
(GHCN) / Climate Anomaly Monitoring System (CAMS) 



Fig. 2. Extrapolation of homogenised De Bilt autumn temperatures 
1901-2005 (crosses) to the value observed in 2006 (horizontal line) 
using a fit to a Gaussian distribution and a Generalised Pareto Distri- 
bution (to the highest 20%). The return times have been computed 
under the obviously false assumption of only interannual variability. 
The upper and lower lines indicate the 95% CI, determined with a 
non-parametric bootstrap. 
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Fig. 3. The temperature anomaly (degrees Celsius) relative to 1961- 
1990) of September-Novemb er 2006 in the GHCN/CAMS dataset 
dFan and van den Doolll2007h . 
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Fig. 4. The return time (years) of September-November 2006 in 
the GHCN/CAMS dataset 1948-2005 under the false assumption of 
a stationary climate, using the more conservative Gaussian extrapo- 
lation. 



Fig. 5. The regression A of local against globally averaged tem- 
perature over 1948-2006. Only grid points where the correlation is 
90% significant (assuming a normal distribution) are shown. 



dFan and van den Dooll 120071) . The GHCN and CAMS time 
series used in the construction of this dataset have not been 
homogenised, in contrast to the De Bilt series of Fig. [Tb 
(which is not included in this dataset). All res ults have been 
verifi ed against the 5° CRUTEM3 dataset dBrohan et all 
20061) (not shown). 

Figure [3] shows that the area with 3°C anomalies stretches 
from the the Alps to Denmark and from Belgium to Poland, 
with maxima in the the Alps and in southern Norway. An 
extrapolation using a Gaussian fit (Fig. |4]l shows that in an 
unchanging climate the return times of this autumn would 
be 10000 years or more in an area shifted somewhat to the 
west of the area with highest anomalies. (The shift is due 
to the smaller variability near the Atlantic Ocean.) As can 
also be shown directly, year-to-year autocorrelations are not 
negligible. 



3 Global warming 

The climate is not stationary: temperatures ha ve been rising 
over Europe as in most of the rest of the world dlPCd.l2007h . 
The effect of this on the probability of the occurrence of an 
anomaly like in autumn 2006 can be studied in a first approx- 
imation by subtracting the local temperature change propor- 
tional to a global temperature, here the 3-yr running mean 
of the global mean temperature, T^^^t). (On interannual 
time scales T^ lohal (t) contains clear ENSO signals, and these 



should not enter in the description of how global temperature 
affects Europe, hence the low-pass filtering.) This gives as a 
first approximation: 



where e(t) denotes the part of the temperature not directly 
proportional to global temperature changes. The coefficients 
A are determined by a fit of local to global temperature up to 
2005 and are shown in Fig. [5] 

Station inhomogeneities (or local temperature anomalies) 
are visible in the GHCN/CAMS dataset (these are not visi- 
ble in the much coarser CRUTEM3 dataset). On average the 
temperature in Europe has increased somewhat faster than 
the globally averaged temperature in autumn, i n accordance 



with the Cold Ocean / Warm Land pattern (ISutton et al. 
2007). For the homogenised De Bilt time series, A = 



T'(t) = AT^ l (t)+e(t), 



(1) 



1.68 ± 0.37 (lcr standard error) over 1901-2005. 

Subtracting the local trend AT^^^t) from the observed 
temperatures, the weather anomalies e(t) remain. The 
anomalies up to 2005 are described reasonably well by a 
Gaussian in autumn; the skewness is —0.46 ± 0.36, so the 
Gaussian fit may even overestimate the positive tail. The au- 
tocorrelation is now slightly negative at lags up to 4 year, so 
that consecutive years of this series can be considered inde- 
pendent. 

In Fig. [6] the cumulative distribution of the autumn tem- 
peratures minus the linear effect of global mean temperature 
changes 1901-2005 is shown. Extrapolating this to the value 
observed in 2006 gives a return time of 650 years. Using a 
GPD to fit the highest 20% of the distribution gives a return 
time of about 3000 years. To compute the statistical uncer- 
tainty margins in these numbers a non-parametric bootstrap 
was performed simultaneously on the regression and the ex- 
trapolation. This gives a 95% CI on the return time of 125 
to 10000 years for the Gaussian extrapolation, 200 years to 
infinity for the GPD extrapolation. The lower bounds are 
reached when the real trend is higher than the central fitted 
value of A = 1.68, lowering the value for 2006 (lower hori- 
zontal line). 
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Fig. 6. As Fig.|2]but now for T(t) - AT^l^i), the De Bilt tern 



perature with the regression against the global mean temperature 
anomalies (with a 3-yr running mean) subtracted. The 95% un- 
certainty interval in the 2006 anomaly above the trend due to the 
uncertainty in the regression is denoted by the horizontal lines. 



The same extrapolation in the GHCN/CAMS dataset 
(Fig. |7]i show the improbable weather to have extended over 
a large part of Europe, with estimates of return times of 
e(t) longer than 200 years over most of the area where the 
anomaly was largest, reaching 2000 years in northern Ger- 
many. The lower bounds of the two-sided 95% confidence 
interval (determined as in Fig. [6]) are more than 300 years 
there. The Gaussian approximation probably underestimates 
the return times, as the anomalies e(t) are negatively skewed 
in most of the area. 

We conclude that global warming has made a temperature 
anomaly like the one observed in autumn 2006 between 10 
and 50 times more likely than under the false assumption of 
stationary climate, with the larger factors near coasts where 
the trend is larger compared to natural variability. Still, other 
factors than global warming conspired to give estimated re- 
turn times well over 200 years in most of the area with large 
anomalies, reaching 2000 years (lower bound 300 years) in 
northern Germany. A rare event indeed, even taking into ac- 
count a shift in the probability density function proportional 
to global temperature changes, determined over the period 
before 2006. 



4 Circulation 

It is obvious that a large part of the temperature anomalies 
in autumn 2006 was caused by the abnormal circulation in 
these months. A persistent low-pressure area over the At- 
lantic Ocean caused predominantly south-easterly winds in 
September, southerly winds in October and south-westerly 
winds in November to the area north of the Alps (Fig. [S). 



In each of these months this corresponded to the direction 
advecting the highest temperatures. 

The question how unusual the circulation patterns were is 
not easy to answer. The local pressure anomalies did not 
reach 3er deviations from the 1961-1990 normals (black lines 
in Fig. [8]). A deviation of this size is not rare without restric- 
tion on the location. 

In order to do a meaningful extreme value analysis on the 
circulation patterns these 2-dimensional fields have to be re- 
duced to local time series. One way to do this is to compute a 
circulation-dependent temperature. In Europe, temperatures 
are determined to a large extend by wind direction due to 
large climatological temperature gradients. Cloudiness also 
affects the temp erature independently. In the very s imple 
model (VSM) of Ivan Ulden and van Oldenborghl d2006l) the 
anomaly in local monthly mean 2m temperature is therefore 
decomposed as 



T'(t) = T c ' irc + Z onchc (t) + MT'(t - 1 

^circ = ^W^westt' 

KonoUt) = A g T^ al (t) +r,(t) 



A s G' south (t) + BV'(t) 



(2) 
(3) 
(4) 



The effect of circulation on the temperature is approximated 
by the circulation-dependent temperature anomalies T c ' irc , 
which are linearly proportional to the zonal and meridional 
geostrophic wind anomalies G' west and G' south and a mea- 
sure for cloudiness, the vorticity V'. The geostrophic wind is 
computed with the corners of a 20° longitude by 10°latitude 
box; the vorticity as the difference between the average of 
the corners and the central value. (All anomalies are relative 
to the mean observed values for 1961-1990.) 

The non-circulation-dependent temperature anomalies 
^noncirc W consist of the part linearly proportional to global 
warming and the remaining noise Finally, M is a mem- 
ory term for past local temperature. This term is modelled 
as a r egression on the previous mo nths' anomalous temper- 
ature (Ivan den Pool and Naplll98ll) . The geostrophic winds 
and vorticity are c omputed from the NCEP/NCAR reanalysis 
sea-level pressure (Kaln ay et al. . 1996 ). 

There is some ambiguity in this model if the cli- 
mate change involves changes in the circulation patterns 
parametrised by the geostrophic wind. However, the in- 
terannual variations in circulation have so far been much 
larger than the long-term shifts, so that in practice the terms 
As, Aw, B are most strongly influenced by the interannual 
variability and this part of climate change is not contained in 
the circulation-independent temperature changes. In autumn, 
there has been no 90% significant change in sea-level pres- 
sures in Europe apart from a slight increase over the Balkan. 

The five model coefficients are fitted at each grid point 
over 1948-2005 (averaged to 1° x 1° for computational rea- 
sons). Averaged over the autumn, the coefficients Aw and 
As reflect the gradients in the climatological temperature 
over Europe (Fig. [9]). In this season the southerly component 
is most important in determining the temperature. Positive 
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Fig. 7. As Fig.|4]but now with the local linear regression against the global mean temperature anomalies (with a 3-yr running mean, Fig. [5} 
subtracted (a), lower bound of the 95% confidence interval on the return time (b). 




Fig. 8. NCEP/NCAR sea-level pressure anomalies (hPa) relative to 1961-1990 in September (a), October (b) and November 2006 (c). The 
black lines indicate the number of standard deviations of the anomaly. 



vorticity leads to more sunshine, which has a positive influ- 
ence on temperature in central and southern Europe, but in 
northern Europe a lack of clouds increases night-time radia- 
tion more and hence cools the surface. The memory term is 
large near seas due to the thermal inertia of sea water on the 
monthly time scale. 

The VSM Eqs (|2]-|4]l explains over half the variance of the 
temperature, i.e., the correlation between the modelled tem- 
perature with rj(t) — and the observed temperature is about 
r = 0.7 to 0.8 in autumn in Europe (FigfTOl). 

In Fig.QT]the contribution of the various terms for the 2006 
event in the VSM are shown. The anomalous circulation con- 
tributed 1.5-2.0°C to the observed anomaly within the linear 
framework of Eqs (|2l|4|. The anomalous vorticity that gave 
rise to large amount of sunshine in September increased the 
temperature by less than 0.5 °C in the seasonal average in this 
linear approximation. 

On the shores of the North Sea persistence has also con- 
tributed. However, this was not due to the below-normal tem- 
peratures in August. The North Sea was still warm from the 



exceptionally high temperatures in July. This is not captured 
by the VSM, which only depends on the previous month, 
hence we cannot give a quantitative estimate. Based on the 
observed SST anomalies of about 2°C at the beginning of 
September it is estimated that this contributed roughly half a 
degree to the autumn temperature anomaly in De Bilt. 

Finally, the return times of T C i rc in autumn 2006 are shown 
in Fig. [12] Comparing this with Fig. [71 the linear effect of the 
circulation in 2006 explains most of the return times in the 
western half of the area with high return times, but in Ger- 
many other effects were responsible for the high anomalies. 
As the memory term is also small in this area, these are not 
described by the simple model. 



5 Climate model simulations 

Observed autumn temperatures are far out of the range ob- 
served so far, even after shifting the distribution with an 
amount proportional to global mean temperature anomalies. 
Do current climate models predict this type of events to hap- 
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Fig. 9. The coefficients of the VSM Eqs j_l4b averaged over September-November: zonal geostrophic wind /IvKpKm s] (a), the meridional 
zonal wind AsfKm. s] (b), the vorticity B[K hPa -1 ] (c); the memory term M [1] in September (d). 
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Fig. 10. Correlation of the September-November temperature 
anomaly from the VSM with rj(t) = and the observed temper- 
ature 1948-2005. 



pen more frequently as Europe heats up? There are two 
caveats when using these models: 

1. the midlatitude circulation response of climate 
models varies greatly from model to model 
and up to now lacks a soun d theo re tical footing 
( van Ulden and van O ldenborgh, 2006; Mille r et al. 



2006); 



2. the local temperature response to global warming is un- 
certain. 

To investigate the changes in the distribution of autumn 
temperatures results from t he 17 standard runs w ith the 
ECHAM5/MPI-OM 1 model djungclaus et all 120061) runs of 
the ESSENCE projecQ were used. These cover the pe- 
riod 1950-2 100 using observed concentrations of greenhouse 
gases and aerosols up to 2000 and the SRES A1B sce- 
nario afterwards. This model simulates the mean circula- 
tion over Europe reasonable well on monthly time scales 



See www.knmi.nl/~sterl/Essence/ for details 
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Fig. 11. Contribution (degrees Celsius) of the geostrophic wind term (a) and vorticity term (b) to the temperature anomaly in autumn 2006. 
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Fig. 12. As Fig. [2] but only for the circulation-dependent tempera- 
ture Tcirc defined by Eq.[3] 



.van Ulden and van O ldenborgh, 2006; van d en Hurk et al. 
2006). In simulations of the 20th century (20c3m experi 



ments) the monthly mean sea-level pressure fields resemble 
those of the observations best of the models in the World Cli- 
mate Research Programme's (WCRP's) Coupled Model In- 
tercomparison Project phase 3 (CMIP3) multi-model dataset. 

An estimate of the systematic errors is provided by a 
comparison with other models that simulate a reasonable 
mean climate over Europe in this measure. These are 



GFDL CM2.1 JPelworth et a 1 120061) . MIRO C 3.2 T106 



(IK- 1 model developersl |2004 . H adGEMl djohnsetal 



2004) and CCCMA CGCM 3.2 T63 dKim et aUl2002l) . The 



MIROC high resolution model did not have enough data to 
reconstruct changes in the full temperature distribution, only 
the mean. 



The ECHAM5/MPI-OM model used in ESSENCE sim- 
ulates the global mean temperature very well, with a ratio 
between observed and modelled trends of 1.10 ± 0.07 (la 
errors). The GFDL CM2. 1 model has similar agreement, but 
the other models overestimate the trend in the global mean 
temperature over 1950-2006 by factors 1.5 (HadGEMl), 1.6 
(MIROC) and 2.0 (CCCMA). To account for these biases, 
we defined the local trend as a regressi on against modelled 



global mean temperature, as was done in Ivan den Hurk et al 
(2006). The local temperature rise as a function of time 
rather than global mean temperature is higher by the same 
factors 1.5 to 2.0 in these models. 

Figure[l3]shows the ratio of observed and modelled warm- 
ing trends in Europe in autumn. The ECHAM5/MPI-OM 
model is seen to underestimate the warming trend in the area 
of the observed extreme by a factor 1.5 or more. In the other 
models the ratio between local observed and modelled tem- 
perature trends has larger errors as there are fewer ensem- 
ble members available, but these models also show a higher 
observed than modelled warming relative to the rest of the 
world in the areas of the autumn 2006 anomaly. 

Figure [14] shows the extreme value distribution of the 
model surface air temperature at the position of De Bilt for 
different 30-year intervals. The distributions are described 
well by a Gaussian. Above the linear increase in tempera- 
ture proportional to the global mean temperature rise, there 
is no indication of any change in the distribution that would 
make extremely warm autumn temperatures more likely. In 
summer this change to more positive skewness is clearly 
seen by steeper slopes in the cumulative distributions (not 
show n); this can be under stand from soil moisture effects 
(e.g., [S char and Jen dritzkv 
Fischer et alll2007t) . 



2004 ISeneviratne et all 120061: 



This result was confirmed for the other climate models 
with a reasonable circulation over Europe for which a com- 
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Fig. 13. Ratio of observed and modelled trends 1950-2006. ESSENCE (ECHAM5/MPI-OM) (a), GFDL CM2.1 (b), MIROC 3.2 T106 
(c), UKMO HadGEMl (d) and CCCMA CGCM3.1 T63 (e). Only grid points where the difference with one is at least one standard error 
(assuming a normal distribution) are shown. The model trends have been computed as a regression against the modelled global mean 
temperature. 
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Fig. 14. Extreme autumn temperatures at 52°N, 5°E in 17 
ECH AM5/MPI-OM 1 20C3M/SRES A1B runs in 1951-1980, 1981- 
2010, 2011-2041, 2041-2070 and 2070-2100. 

parison between the 22nd and 23rd century with the 20th 
century could be made. None of these show an increase in 
the slope of the cumulative distribution at the grid point cor- 
responding to De Bilt, Fig. IT31 

The probability of temperatures as high as those observed 
in 2006 increases if global warming shifts the distribution 
to higher values, or the width of the distribution increases 
towards higher temperatures. The climate models considered 
here show no evidence for either of these mechanisms: they 
simulate a lower shift of the distribution and no change in 
shape. Note that the change in width is impossible to obtain 



from observations. 
6 Conclusions 

The autumn of 2006 was extraordinarily warm in large parts 
of Europe, with temperatures up to 4°C above the 1961-1990 
normals. Assuming an unchanging climate, this would cor- 
respond to return times of 10000 years and more. 

Global warming has made a warm autumn like the one ob- 
served in 2006 much more likely by shifting the temperature 
distribution to higher values. Taking only this mean warm- 
ing into account, the best estimate of the return time of the 
observed temperatures in 2006 still is more than 200 years in 
large parts of Europe. The lower bound is 100 years or more 
years, reached when the trend is much larger than estimated 
from the years before 2006. 

Apart from global warming, the anomalously high tem- 
peratures in Europe in autumn 2006 were mainly caused by 
the linear effects of a persistent southerly wind direction ad- 
verting warm air to the north and persistence from the very 
hot July along the shores of the North Sea. A simple lin- 
ear model incorporating these factors reproduces most of the 
temperature anomalies west of the Alps, but not in Germany. 

Current climate models already underestimate the ob- 
served mean warming in Europe relative to global warming 
before 2006. They also do not show an extra increase of the 
warm tail of the distribution as the climate warms. Either au- 
tumn 2006 was a very rare event, or these climate models do 
not give the correct change in temperature distribution as the 
temperature rises. 
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